Mon. Not. R. Astron. Soc. 000, ITHT51 (2010) Printed 7 June 2010 (MN LJO^X style file v2.2) 



A Large Sample of Photometric Rotation Periods for FGK 
Pleiades Stars 



O ; J. D. Hartman 1 *, G. A. Bakos 1 , G. Kovacs 2 , R. W. Noyes 1 

' _ ^ , 1 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 
O , > KmMy Observatory, Budapest, Hungary. 

7 June 2010 



6 



> 
o 
in 

On 
O 

\6 
O 
O 



X 



ABSTRACT 

Using data from the HATNet survey for transiting exoplanets we measure photomet- 
ric rotation periods for 368 Pleiades stars with 0.4 Mq < M < 1.3 Mq. We detect 
periodic variability for 74% of the cluster members in this mass range that are within 
our field-of-view, and 93% of the members with 0.7 M Q < M < 1.0 M Q . This in- 
creases, by a factor of five, the number of Pleiades members with measured periods. 
Included in our sample are 14 newly identified probable cluster members which have 
proper motions, photometry, and rotation periods consistent with membership. We 
compare this data to the rich sample of spectroscopically determined projected equa- 
torial rotation velocities (vsmi) available in the literature for this cluster. For stars 
with M > 0.85 Mq the rotation periods, vsini and radius estimates are consistent 
with the stars having an isotropic distribution of rotation axes, if a moderate differen- 
tial rotation law is assumed. For stars with M < 0.85 Mq the inferred sini values are 
systematically larger than 1.0. These observations imply that the combination of mea- 
sured parameters P(v sin i)/R is too large by ~ 24% for low- mass stars in this cluster. 
By comparing our new mass-period relation for the Pleiades to the slightly older clus- 
ter M35, we confirm previous indications that the spin-down stalls at ~ 100 Myr for 
the slowest rotating stars with 0.7 Mq < M < 1.1 Mq - a fact which may indicate 
that the internal transport of angular momentum is inefficient in slowly rotating solar 
mass stars. 

Key words: stars: rotation, spots, late-type — open clusters and associations: indi- 
vidual: Pleiades — techniques: photometric catalogues 



1 INTRODUCTION 

Due to its proximity, richn ess, and age, the Pleia des star 
cluster (d ~ 133 pc, e.g. Soderblom et al.l 120051 ; age ~ 
125 Myr, e.g. lStauffer et al.|[l998l ) has served for many years 
as one of the benchmark cl usters in studies o f stellar an- 
gular momentum evolution. ISkumanichl (|l972n first deter- 
mined the empirical spin-down relation for solar-type stars 
(a; oc t -1 / 2 ) by comparing measurements of the projected ro- 
tation velocities {vsmi) of solar-mass stars in the Pleiades 
and the Hyades to the equatorial rotation velocity of the 
Sun. Since then there have been numerous studies presenting 
rotation periods and vsmi values for stars in the Pleiades 
(see the references given in sections 14. II and 14. 2\ . 

Recently, large samples of stellar photometric rotation 
periods have been published for a number of open clus- 
ters (see for example the summary and references given by 
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llrwin fc Bouvied 120091 for 18 open clusters, together with 
some new re sults, not included in their list, for Coma 
Be renices by ICollier Cameron et all I2009I . for NGC 2301 
by ISukhbold fc Howelll 120091 . and for M34 by Ijames et all 
2010). These data have enabled theoretical investigations of 
the angular momen tum evolution of low-m ass stars (again 
see the review by llrwin fc Bouyierl |2009| . and also a re- 
cent study by iDenissenkov et aL 20091 ). which have led to 
a number of insights into stellar physics. Observations of 
very young clusters with ages < 10 Myr can be used to 
study phenomena suc h as ac cretion-driven stellar winds 
dMatt fc Pudrit d 120051 l2008al lbh or star-disk-locking (e.g. 
lKonigllll99ll : IShu et al.lll994l ). Observations of clusters with 
ages ~ 50 — 200 Myr are essential for constraining the 
time-scale of coup ling between s tellar convective and ra - 
diative zones (e.g. iBouvierllioog ; IDenissenkov et al.l 120091 ) . 
Finally, observations of older clusters are imp ortant for con - 
straining models of magnetized stellar winds ( Kawalcr 1988) 
and for calibrating the rotation-age relation (e.g. Barnes! 
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Figure 1. Left: A typical 10.6° X 10.6° HAT-9 image of HATNet field G259 containing the Pleiade s cluster. The 2° and 5° radius circles 
contain 68% and 100% of the candidate cluster members in the catalogue of lStauffer et al~l | |2007| ). Right: A 1.'5 X 1'5 zoomed-in view 
of the center of the cluster. 



20071 : iMamaiek fc Hillenbrandll200Sl : ICollier Cameron et al.l 
20091 ). 



Although the comparison of models to the observed ro- 
tation period distributions in open clusters has led to a num- 
ber of insights, these studies are hindered by incompleteness 
and poorly determined selection effects in the available data- 
sets. These effects may in turn lead to incorrect theoretical 
inferences. The available photometric rotation periods for 
the Pleiades, in parti cular, show evidence of being biased 
toward short periods (|Denissenkov et al]|2009l ). this is un- 
surprising since longer period stars generally have lower am- 
plitude variations and also require observations spanning a 
longer time-base to be detected, moreover apparent varia- 
tions due to systematic errors in the photometry may domi- 
nate on long time-scales; the results to be discussed here do 
not suffer from this bias because of the long time-base and 
high precision of our observations. It is only possible to de- 
termine whether or not the photometric periods are biased 
for this cluster because it is fairly unique in having a rich, 
unbiased sample of usini values in the literature. Given the 
importance of the Pleiades for the study of stellar angular 
momentum evolution, an unbiased sample of rotation peri- 
ods for stars in this cluster would be quite useful. 

In this paper we present photometric periods that we 
associate with rotation periods for 368 Pleiades stars. This 
increases the number of Pleiades stars with measured peri- 
ods by a factor of 5, and importantly, it is 93% complete for 
stars in the mass range 0.7 M© < M < 1.0 Mq. For this 
survey we use data from the Hungarian-made Automated 
Telescope Network (HATNet fBakos et af1l2004 ), a project 
that uses a network of 6 small-aperture, wide-field, fully- 
automated telescopes to search for trans iting exoplanets or- 
biting bright stars fe.g. lBakos et al . 2010). A similar study of 
stellar rotation for field K and M dwarf stars using HATNet 
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Figure 2. The number of observations of HATNet field G259 
obtained on each night with the HAT-9 and HAT-10 instruments. 



has be en recently presented by our group (Har tman et al.l 
l2009bl ). 

The structure of the paper is as follows: in section [5] 
we briefly describe the photometric observations and data 
reduction; in section [3] we select periodic variable stars, es- 
timate the errors on the period, correct the light curves for 
distortions, present the catalogue of rotation periods, and 
select new cluster members; in section [4] we compare our 
data to previous rotation period and vsini measurements, 
and we also compare our results to other open clusters. We 
summarize the results in section [5] 
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Figure 3. Unbiased RMS (eq. [2} vs. Sloan r magnitude for the 
EPD/TFA corrected light curves of stars in field G259 (gray-scale 
points). The dark filled circles show probable Pleiades members 
(section I3.4|l . Stars brighter than r < 10 are saturated in a signif- 
icant fraction of the images. 



2 OBSERVATIONAL DATA AND INITIAL 
REDUCTION 

Photometric time-series observations of the Pleiades cluster 
were obtained between 15 September 2008 and 16 March 
2009 using the identical HAT-9 and HAT- 10 11 cm aper- 
ture robotic telescopes located at Mauna Kea Observa- 
tory (MKO) in Hawaii and at Fred L. Whipple Observa- 
tory (FLWO) in Arizona respectively. Each telescope used a 
4Kx4K CCD and a Sloan r filter to observe a 10.6° x 10.6° 
field of view (FOV) centered at 03:30:00, 22:30:00 (J2000). 
The Pleiades are not centered in this field; the field (in- 
ternally designated as G259) was observed as part of the 
standard operations of the HATNet transit survey, and was 
not specifically chosen to observe the cluster. Figure [T] shows 
a typical image of this field. The images have a pixel scale of 
9"; the full- width at half-maximum (FWHM) of the point- 
spread function (PSF) is ~ 2.6 pixels (23")- A total of 3648 
exposures of 5 minutes, taken at 5.5 minute cadence, were 
obtained with the HAT-9 instrument, while 3138 images 
with the same exposure time and cadence were obtained 
with the HAT-10 instrument. Figure [2] shows the tempo- 
ral distribution of observations obtained with each of these 
instruments. 

The images were calibrated and reduced to light curves 
using too l s developed for th e HATNet transit survey (see 
[Pal 120091 ; feakos et all |2004 ). Briefly, after applying bias, 
dark current, and twilight sky-flat calibrations in a stan- 
dard fashion, stars were identified on the images. The star 
lists were the n matched to the Tw o Micron All-Sky Sur- 
vey (2MASS; ISkrutskie et al]|2006l ) po int source catalogue 
(PSC) using the methods described bv lPal fc Bakosl (| 20061 ) 
to determine the astrometric solutions for the images. Aper- 
ture photometry is then performed at the positions of all 
2MASS sources with r < 14.5 transformed to the image 
coordinate system. For each resulting light curve the me- 
dian magnitude is fixed to the approximate r magnitude 
of the source based on the following tra nsformation from 
the 2MASS J, H and K s magnitudes (see lBilir et~al"1l2008l . 
who give transformations from 2MASS magnitudes to Sloan 



(g — r) and (r — i) colors; we followed the same procedure 
to determine a transformation from 2MASS magnitudes to 
Sloan r magnitude): 



0.6975 + 2.9782J - 0.8809// - 1.1230/sTs. 



(1) 



The RMS scatter between the observed and predicted r 
magnitudes for the stars used to determine this relation is 
0.076 mag. For each source, photometry is performed using 
three apertures of radii 1.45, 1.95 and 2.35 pixels. Follow- 
ing the post-processing routines discussed below, we adopt 
a single "best" aperture for each light curve. For stars with 
r > 13.5 we take the smallest aperture to minimize sky 
noise, for brighter stars we adopt the aperture for which the 
processed light curve has the smallest RMS. 

The initial ensemble calibration of the light curves 
against variations in the flux scale w as perform ed using the 
method described in section 2.7.3 of (Pal 2009j). These light 
curves are then passed through two routines that filter out 
systematic variations that are correlated with measureable 
instrumental parameters or are present in other stars in 
the field. The first routine, external parameter decorrela- 
tion (EPD), decorrelates each light curve against a set of 
measured instrumental parameters including the shape of 
the PSF, the sub-pixel position of the star on the image, 
the zenith angle, the hour angle, the loc al sky background , 
and the variance of the background (see lBakos et al.ll2010l ). 
This decorrelation is done independently on the data from 
the HAT-9 and HAT-10 telescopes. The procedure is applied 
assuming that each star has a constant magnitude. For large 
amplitude variable stars this may distort the signal and may 
lower the S/N of the detection, however such large ampli- 
tude variables will generally still be detectable. 

After applying EPD, the light curves are then 
processed with the Trend -Filtering Algorithm (TFA; 
iKovacs. Bakos fc Novesll2005l ) which decorrelates each light 
curve against a representative sample of other light curves 
from the field. We used 530 template stars (~ 8% of the 
total number of images for the field). In applying the TFA 
routine we also clip 5a outliers from the light curves. At 
this point in the analysis we apply EPD and TFA in signal- 
recovery mode (i.e. we apply them under the assumption 
that the signal is constant) , once a signal is detected we then 
apply EPD and TFA in signal-reconstruction mode on the 
original light curve to obtain an undistorted trend-filtered 
light curve for the star (see section l3~3"f) . As for EPD, signal- 
recovery mode TFA may distort the signal and lower the 
S/N of large amplitude variable stars, though typically the 
effect is not significant enough to prevent the variable from 
being selected. 

In figure[3]we show the unbiased RMS of the EPD/TFA 
corrected light curves for all stars in field G259. The unbi- 
ased RMS of a light curve is calculated using 



RMS 



Mr 



N-Np 



(2) 



where rm are the individual magnitudes, (m) is the av- 
erage magnitude, N is the number of points in the light 
curve, and Np — 544 is the number of parameters used 
in applying EPD/TFA. We mark separately the probable 
Pleiades members (see section T3.4p . Stars with r < 10 are 
saturated in a significant fraction of the images (due to vi- 
gnetting the exact magnitude of saturation depends on the 
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position of the star on the image; it also depends on the 
sky brightness and transparency which changes from im- 
age to image). The unbiased RMS of stars with r ~ 10 is 
~ 4 mmag. The light curves of Pleiades members generally 
have greater RMS than those of field stars; this reflects the 
fact that detectable photometric variability is significantly 
more common for Pleiades members than it is for field stars 
(see section [3~4l) , 



3 SEARCH FOR PERIODIC VARIABLES 



We use the Lomb-scargle pe riodogram l|Lomblll976l : IScargld 
1 19821 : IPress fc RvbickJll989h as implem ented in the VAR- 
TOOLS program (|Hartman et alj|2008bl ) to search the light 
curves for periodic variations. We generate the periodogram 
for each light curve at a frequency resolution of 0.1 /T be- 
tween 0.01 d _1 and 10.0 d —1 , where T is the time-span of a 
given light curve. The high-frequency cut-off of 10.0 d -1 is 
adopted as stars are not expected to have rotation periods 
shorter than this, and expanding the range of the frequency 
search increases the bandwidth penality, decreasing the sig- 
nificance of a given detection. We correct the periodogram 
for red noise by removing any low-order frequency depen- 
dence in the mean value of the periodogram and in the RMS 
of the periodogram. To do this we fit a fifth order polyno- 
mial to the periodogram, as well as a fifth order polynomial 
to the RMS of the periodogram calculated in 100 frequency 
bins. We then define a new periodogram using 



p'(w) 



(p(w)-p(o;)) + (p> 



(3) 



where p is the value of the original periodogram at frequency 
w, ct(uj) is the polynomial fit to the RMS of the periodogram 
as a function of frequency, p(cj) is the polynomial fit to the 
periodogram, and the brackets denote averaging over fre- 
quency. We perform the fit using 5-cr iterative clipping, using 
the model value of <r(w) in the clipping. 

We identify the highest peak in the corrected peri- 
odogram p' and determine its S/N using an iterative 5-cr 
clipping in calculating the RMS of the periodogram. To set 
the selection threshold we simulate 10000 white noise light 
curves with the same time sampling as the observations, 
calculate the red-noise corrected L-S periodogram of each 
simulation, and determine the S/N of the highest peak in 
each periodogram using the same procedure as used for the 
real light curves. Figure [4] shows the false alarm probability 
(FAP) as a function of S/N determined from the simula- 
tions, which we find to be well-fit by a function of the form: 



FAP = 1 - (1 - , 



-0.84(S/N+1K 1240 



(4) 



This is similar to the form expected for a normalized L- 
S periodogram w ithout red-noise correction or clipping (e.g. 
IPress et al.ll 19921 '). We adopt a cut-off of S/N> 16 which cor- 
responds to a 0.08% FAP. A total of 2236 of the 36,011 stars 
in the field pass the selection. In section [3741 below we focus 
on the stars that are likely to be cluster members. Adopt- 
ing a lower threshold of S/N> 13 (corresponding to a 1% 
FAP) would increase the total number of selected potential 
variable stars to 2878. 
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Figure 4. False alarm probability vs. S/N of the peak in the 
red-noise corrected L-S periodogram. The solid line shows the 
results from conducting simulations of white noise light curves, 
the dashed lines shows eg. l4l 



3.1 Resolving Aliases 

In some cases there is an ambiguity in choosing the correct 
period from among several possible aliases or harmonics. In 
general we follow the convention of adopting the highest 
S/N peak in the red- noise corrected L-S periodogram. We 
found, however, that often the red-noise correction changes 
the most signficant peak in the periodogram to a high fre- 
quency alias | f s — fo of the most significant peak in the un- 
corrected periodogram, fo- Here / s is the sidereal frequency. 
In these cases we adopt the period from the uncorrected pe- 
riodogram. We found that generally the periods from the 
red-noised corrected periodograms pile-up near 1 day while 
the periods from the uncorrected periodograms fall along 
the main period- mass relation (see figure [T2J. 



3.2 Estimation of Period Uncertainties 



As discussed by lHartman et al.l (|2009al ) , the primary fac- 
tors that contribute to errors in the measured photometric 
rotation period are: 



Instrumental Effects 

(i) choosing an alias or harmonic of the true period, 

(ii) noise in the photometry together with finite sampling 
of the light curve, 

(iii) inadequacies in the model used to determine the pe- 
riod (e.g., the light curve is periodic but not sinusoidal), 



Physical Effects 

(iv) spot evolution, and 

(v) differential rotation. 

The periodogram for a light curve will contain several 
discrete peaks at harmonics and aliases of the true period. 
Deciding which peak is the correct one is difficult to do in 
general. We do not consider this effect in estimating the 
uncertainties on the measured periods. 
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The effects (ii)|(iv) listed above all contribute to the 
spread in the periodogram peak, which reduces the preci- 
sion with which the position of the peak may be measured. 
We determine the resulting error in the period by fitting 
a Gaussian function to the points within 2/T of the peri- 
odogram peak and measuring the lcr spread in the position 
of the peak us ing the Ex tended Markov-Chain Monte Carlo 
technique fsee |Palll2009f l. We set the error for each point in 
the periodogram equal to the RMS of the full periodogram 
determined with an iterative 5a clipping. We also record the 
standard deviation of the best-fit Gaussian, as this contains 
information on the spot life-times. 

Differential rotation may contribute to the uncertainty 
in two ways. If there are multiple spot groups at different 
stellar latitudes, the differing rotation periods for each group 
will broaden the periodogram peak, and will contribute to 
the uncertainty in the period as just described. Differential 
rotation will also lead to a systematic error in determining 
the rotation period of the star since the degree of differential 
rotation and the latitude of the dominant spot group are not 
known, making it impossible to relate the measured period 
to the period at a reference latitude (such as the equator). 
If Pleiades stars exhibit solar-like differential rotation (the 
rotation period increases towards the poles), the measured 
periods will be systematically longer than the equatorial pe- 
riod. The Sun exhibits differential rotation following 



(*)• 



(6) 



Pp = P E Q/(l-fcsin 2 P) 



(•») 



where Pp is the rotation period at latitude /3, Peq is the 
equatorial rotation p eriod, and k = 0.1 9 when using spots 
to track the rotation (|Kitchatinovll2005l ). For younger, more 
rapidly rotating sta rs, the value of k is expected to decrease 
|Brown et alj|20u4 ). Observations of the P = 8.77 d solar- 
like star k 1 Ceti by the Microvariability and Oscillations of 
Stars (MOST) sat ellite bear this theor etical expectation out, 
finding k = 0.09 (| Walker et alj|2007h . The observations of 
k Ceti also indicate that spots on rapidly rotating stars may 
be found at any latitude; this is in contrast to the Sun where 
spots are rarely seen with |/3| > 30°. Assuming the dom- 
inant spot groups are isotropically distributed on Pleiades 
stars (i.e. sin /3 is uniformly distributed), and that these stars 
exhibit solar-like differential rotation with k = 0.09, we ex- 
pect the measured rotation periods to have a mean value of 
1.03Peq, with a standard deviation of 0.03Peq- For k — 0.19 
the respective values are 1.07Peq and 0.07Peq- We do not 
include systematic errors due to differential rotation in the 
period uncertainties reported in our catalogue. 



3.3 Signal Reconstruction with EPD/TFA 

Once a period is determined for a star, we obtain a new 
trend-filtered light curve by running EPD and TFA in sig - 
nal reconstruction mode fe.g. lKovacs. Bakos fc Novesl2 005). 
This correction is important to get an unbiased measure- 
ment of the amplitude of photometric variations. To do this, 
we fit to the pre-EPD light curve (i.e. the light curve has 
been corrected for ensemble variations in the flux scale, but 
has not had any other filtering applied to it), a model of the 
form: 

10 

m(t) — ap + S (ajCos(2TYit/P) + bjsin(2nit/ P)) 



The first sum is a Fourier series with period P (the period 
found for the light curve in section [3} which is used to repre- 
sent the physical signal, while the second sum is a model for 
the instrumental/atmospheric variations. The free parame- 
ters in this model are the ctj, bi and a coefficients, while the 
Si(t) terms, which represent instrumental variations, con- 
sist of 530 template light curves, together with 14 known 
instrumental parameter sequences (including: the X and Y 
subpixel positions of the star to first and second order, 3 
parameters describing the PSF shape to first and second or- 
der, the hour angle, the zenith angle, the background, and 
the deviation of the background). Because all free parame- 
ters enter linearly in eq. |S] the fit can be done quick ly using 
singular value decomposition (e.g. IPress et al.lll992T ). 

We take the amplitude of the light curve to be equal to 
the peak-to-peak amplitude of the Fourier series in eq. [6] 
For four stars, HAT-259-0000690 (PELS 020), HAT-259- 
0000923 (HII 120), HAT-260-0000924 (HII 1182) and HAT- 
260-0007928 (HII 906), with periods very close to 3, 4, 3 and 
2 days respectively, the poor phase coverage of the HAT- 
Net light curves means that the high order harmonics in 
the Fourier series are not well constrained. For these stars 
we therefore do not include any harmonics other than the 
fundamental (i.e. the Fourier series includes one term only) 
when fitting equation [6] 

The median value of the ratio of amplitudes from signal 
reconstruction EPD/TFA light curves to amplitudes from 
signal recovery EPD/TFA light curves is 1.3. This illustrates 
the importance of using signal reconstruction EPD/TFA to 
get an unbiased amplitude measurement. 



3.4 Catalogue of Rotation Periods for Pleiades 
Members 

Using the recent compilatio n of probable Pleiades mem- 
bers bv lStauffer et al. I l|2007l . hereafter S07), we identify 470 
probable members with magnitudes in the range 9.5 < r < 
14.5 for which we have obtained light curves. We detect pe- 
riodic variations for 350 of these stars (i.e. 74% of them); 
for stars with 11 < r < 13, our period detection rate is 
93%. Faint stars with 13 < r < 14.5 have low photometric 
precision making it more difficult to detect their variations, 
while hotter stars with 9.5 < r < 11 have lower amplitudes 
of variation. Stars with 9.5 < r < 10 may also be saturated 
in a significant fraction of the images. As a comparison, for 
stars with 11 < r < 13 that are not selected as probable 
cluster members, the period detection rate is 9.8%. Two fac- 
tors that lead to a much higher incidence of variability for 
Pleiades stars include their young age relative to field stars 
(so that they are more active and have higher amplitudes of 
variability), and the fact that the majority of field stars are 
more distant than the Pleiades, so Pleiades stars of a given 
mass are brighter, and have higher precision light curves, 
than most field stars of that mass in the survey. 

Our final catalogue of rotation periods is included in 
the supplementary material to the online edition of this ar- 
ticle. Table [T] shows a portion of the catalogue for guidance 
regarding its form and content. We include in the catalogue 
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Figure 5. Phased light curves for 12 stars randomly selected from 
our catalogue, sorted by rotation period. The gray-scale points 
show all the photometric data, the dark filled circles show the 
phase-binned l i ght cu rve. The name for each star is taken from 
Stauf fer et al. I (120071 ) and the listed period is in days. These light 
curves have had signal-reconstruction EPD/TFA applied (Sec- 
tion E3J. 



fifteen stars with rotation periods in the literature that we 
do not recover, including the star CFHT-PL 8 which was not 
included in the S07 catalogue, and exclude one star with a 
period in the literature that is not a cluster member (see 
Section 14.11 for further details on these stars). We also in- 
clude 18 variable stars which are not in the S07 catalogue, 
but which have photometry, proper motions, and periods 
consistent with cluster membership (see section [331 below) . 

Figure [5] shows phased light curves for 12 stars ran- 
domly selected from our catalogue, while figure [5] shows L- 
S periodograms for 3 of these stars. Figure [6] also shows 
an example of a star which is not selected after applying 
the red-noise correction to the periodogram (equation [3]). 
The photometric light curves for all 488 probable mem- 
bers observed with HATNet will be made publicly avail- 
able from the NASA S tar and Exoplanet Database (NStED; 
Ivon Braun et al.ll2009l fl. 

In addition to the catalogue of rotation periods, we also 
provide a catalogue of Pleiades stars observed with HATNet 



http: / /nsted. ipac.caltech.edu 
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Figure 6. L-S periodograms for 3 of the stars displayed in fig- 
ure \E\ and for a fourth field star which passes the S/N selec- 
tion before applying the red-noise correction to the periodogram 
(equation [3]l , but does not pass the selection after applying the 
correction. For each star we show the original L-S periodogram 
on the left, and the red-noise corrected periodogram on the right. 
We list the period (in days) and the frequency (in cycles per day) 
selected for each star. For the star in the bottom panel, this is 
the period that is selected from the un-corrected periodogram. 
For HCG 353 the low frequency alias is adopted following the 
procedure described in Section |3~T1 



for which we do not detect a period (Table [2} , a catalogue 
of field variable stars (Tabled, and a list of previous period 
measurments for Pleiades stars (Table [4]). 



3.5 New Cluster Members 

Here we leverage the enhanced photometric variability of 
Pleiades stars relative to field stars to identify new mem- 
bers of the cluster. A similar method for select ing Pleiades 
members was employed by lHaro et all l|l982l ) who iden- 
tified flare stars as candid ate members. More recently, 
ICollier Cameron et al ] (|2009l ) used photometric variability 
to select members of the Coma Berenices open cluster. 

Figure[7]shows the selection of cluster members. We first 
select variable stars that are within 5° of the centre of the 
cluster, and are near the cluster main sequence on a J — Ks 
colour-magnitude diagram (CMD), taking the photometry 
from the 2MASS 6x PSC (see the description by Cutri et 
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alQ) when available, and from the main 2MASS PSC for 
all other stars. A total of 221 of the 958 variables not in- 
cluded in the S07 catalogue that are within 5° of the cluster 
centre pass this selection. Of these, 18 9 have proper mo - 
tions from the PPM- Extended catalogue (|Roser et alj [2008). 
while 28 of the remaining 3 2 stars have proper motions from 
the USNO-B1.0 catalogue (|Monet et aJj|2003h . To establish 
the proper motion membership probabi l ity w e follow the 
procedure outlined in lDeacon fc Hamblvl (|2004T ). Briefly, we 
assume that proper motions are distributed as: 



* = /*/ + (!-/)#_ 



(7) 



where / is the fraction of stars in the field, <£>/ is the distri- 
bution for field stars, and $ c is the distribution for cluster 
members. We assume bivariate Cauchy distributions for the 
field and the cluster of the form: 



2*rj (i + ((^ - f, xf )2 + _ M f)/vf) 

__i i 

2-nVl (i + ((^ - Ma;c )2 + ( Ms/ - /v0 2 )/r§) 



(8) 
(9) 



where fj, x and /j, y are the proper motion in the right ascen- 
scion and declination directions respectively. We determine 
the free parameters /, fj, x f, fi v f, (i xc , fi yc , Tf and P c by using 
the downhill simplex algorithm to maximize the likelihood 
function 



L = ^ln$, 



(10) 



where the sum is over all variable stars. We also attempted to 
use Gaussian and Exponential distributions for <&j> and $ c , 
but found that the Cauchy distributions provided the best 
match to the observations. The membership probability for 
a star is then given by 



P = 



(1 - /)*c 



f* f + (1 - /)*„ 



(11) 



We require p > 70% for the star to be considered a probable 
member. Of the 22 stars which pass the above selections, we 
reject 4 with periods longer than any known cluster members 
of comparable magnitude. 

Four of the 18 remaining probable cluster members 
have previously been selected as potential cluster members, 
but were not included in the S07 catalogue. These include 
HCG 84 and HCG 235, two well-known flare stars include d 
in the catalogue of Pleiades members bv lHaro et al. 

I (|1982| ). 

and in several subsequent membership catalogues, and the 
stars SRS 79807 and SRS 34337 which we r e iden tified as 
proper-motion members by ISchilbach et~al1 (|l995T ). but to 
our knowledge, have not appeared in subsequent studies. 
All 18 stars are included in the rotation period catalogue. 



4 DISCUSSION 

4.1 Comparison With Other Period 
Measurements 

To compare our measurements to previous observations of 
the Pleiades, we use the compilation of rotation periods 
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Figure 7. Top: J — Kg CMD for all variable stars. Cluster 
members from S07 are plotted with filled circles, variables with 
J — Kg photometry consistent with cluster membership are plot- 
ted with open circles, and variables that are not photometric 
members are plotted with crosses. The solid lines dilineate the 
selection of candidate cluster members. Centre: Vector Point Di- 
agram (VPD) for all identified variable stars. Filled circles show 
cluster members from S07, crosses are stars which were rejected as 
non-photometric members in the left panel, open circles show pho- 
tometric members which also have a proper motion membership 
probability > 70%, Xs show photometric members which have a 
proper motion membership probability < 70%. Bottom: Period 
vs. ify-magnitude. Filled circles show cluster members from S07, 
open circles correspond to new candidate members which passed 
the other two selections and also have periods that are consistent 
with cluster membership. Crosses show candidate members that 
passed the previous selections, but have periods that are too long 
to be cluster members. 



2 http: / /www. ipac.caltech.edu/2mass/releases/allsky/doc/seca3_l. html 
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provided on the WEBDA databas^fl The original measure- 



ments come f r om a v arie ty of so urces dStauffer et al. Ill987l: 



Prosser et all Il993allbl Il995l: IMarilli. Catalano fc Frascal 
1997; Krishn amurthi et al 
Messinal l200ll; 



1998; Tcrndrup ct al 



C lar ke. MacDonald. fc Owen; 



199£ 
2004 



Scholz fc Eisloffell 120041). We also include period s for 
11 stars from van Leeuwen. Alphenaar. fc Meva (jl987) that 



are not listed on WEBDA. A total of 66 stars are included 
in the compilation, which is provided in tabled] 

We recover rotation periods for 50 of the stars with 
previous measurements. Here we discuss the 16 stars 
with previous measurements which we do not recover. 
The 9 very low-mass stars studied by IScholz fc Eisloffel 
l|2004l ) are all too faint to be observed with HATNet, 
the st ar HD 23386 st u died b y IMarilli. Catalano fc Frascal 
i|l997h and [m 

cssm 2 (|200ll ) is saturated in our im- 
ages. There are four stars that we observed but 

did not det ect a period for, including: H II 708 

ilProsser et al.l 1 1993bl: IMarilli. Catalano fc Frascal Il997l) 



HII 727 (iProsser et all 1 1993bl: IMarilli. Catalano fc Frascal 
Il997h . HI I 975 (IMarilli Catalano fc Frascal Il997l ) and 
HHJ 409 (|Terndrupet al.l ll999l ). The first three of these 
stars are near the bright end of the magnitude range 
covered by HATNet, while HHJ 409 is near the faint 
end. Finally there are two stars listed as variables on 
WEBDA that are not included in the S07 catalogue. One 
of these stars, HII 3167, was r eport ed to be variable 
by IClarke. MacDonald. fc Owens! i|2004h with a 0.25 mag 
amplitude in B, 0.15 mag amplitude in V and a pe- 
riod of ~ 0.9 d. We find no evidence for variabil- 
ity in our light curve for this star (we match the star 
to 2MASS J2000.0351555+2442326 based on the plate 
x/y position given on WEBDA, however this identifica- 
tion may not correspond to the s ame star observed by 
IClarke. MacDonald. fc Owens! 12001 ) . In any case the star 
is not a cluster member (based on the color a nd magnitude 
given bv lClarke. MacDonald. fc Owensll2004l) . so we do not 
include it in the final catalogue of Pleiades variables. The 
other known variable that is no t in the S07 catalogue is 
CFHT-PL 8 dBouvier et al.ll 19981) . which was identified as a 
variable by iTerndrup et"a!T(ll999l ). This very low mass star 
is too faint to be observed by HATNet. We include it in 
the final catalogue of rotation periods, with coordinates and 
Ks photometry taken from the 2MASS catalogue, and V 
magnitude taken from lSamus et al. I (|2003h . 

For 44 of the 50 previously studied stars that we re- 
cover, our frequency is within 0.1 d -1 of at least one of 
the previously published frequencies for the star; and for 35 
stars our frequency is within 0.01 d _1 of at least one of the 
previously published frequencies. Figure [8] shows a compar- 
ison of the rotation frequencies measured in our survey to 
those given in the literature for stars with previous measure- 
ments. Below we discuss the 6 cases where our frequencies 
are substantially different from all of the literature values: 

HII 250 - Has publis h ed pe riods of 0.843 days from 

Marilli, Catal ano fc Frascal (|l997l ) and 0.591 days from 
Messinal i|200ll ). We measure a period of ~ 4.232 days. Our 



light curve shows no evidence for periodicity at either of the 
published periods. 



http: / /www. univie.ac.at/webda/webda. html 
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Figure 8. Comparison between the rotation frequencies mea- 
sured in our survey and those given in the literature for stars 
with previous measurements. For stars with multiple frequencies 
given in the literature, we show the literature value that is closest 
to our own (filled circles); the open circles show the other liter- 
ature values for these stars. Points between the solid lines have 
literature frequencies that are within 10% of our frequencies. The 
dotted lines show the ±1 d~ 1 alias frequencies, and the 2.0 and 
0.5 harmonic frequencies. 



335 



HII 

from Stauffcr ct al 



Ha s pub lished periods of 0.2 65 days 
dl987h a nd 0.4073 days from IMessinal 



(2001). Stauffcr ct al. (1987) also report an alias period of 
0.360 days for this star that is consistent with our measured 
period of 0.357 days. 

HII 885 - The only published period for t his star is 
0.435 days from lMarilli. Catalano fc Frascal (|l997l ). We mea- 
sure a period of 6.83 days, and find no evidence for period- 
icity at 0.435 days. 

HII 1039 - The only publis hed period for this star is 
0.784 days from IMessinal |200ll ). We measure a period of 
5.22 days, and find no evidence for periodicity at 0.784 days. 

HI I 1124 - Ha s publi shed periods of 5.9 days 
from Prosser et al. (Il995l) a nd 6 days from 
Ivan Leeuwen. Alphenaar. fc Mevsl (1 19871 ). We find a 
period of 0.858 days for this star, however prior to cor- 
recting the periodogram for red-noise the top peak in the 
periodogram for this star is at 6.133 d ays. This star has 
wsini = 3.50 kms -1 (|Queloz et al.|[l99^ ). which is closer to 
the expected equatorial rotation velocity of ~ 5.8 kms - 
for the ~ 6 day period than the velocity of ~ 40 kms -1 for 
the shorter period. We therefore adopt the longer period 
for the final catalogue. 

HII 1653 - Both iKrishnamurthi et all (|l998h and 
IMessinal (|200ll ) measured periods of ~ 0.75 days for this star. 
Our period of 0.368 days is half the value of the published 
periods, and is inconsistent with the previously published 
light curves for the star. We therefore double our period for 
the final catalogue. 



4.2 Comparison with vsmi data 

We compare our period measurements to the rich sample of 
vsmi measurements available for members of the Pleiades. 
Measu r ements of v sin i are t aken from IStauffer et al.l 
l|l987h , ISoderblom et~aH (|l993l ). lOueloz et all (|l998l ). and 
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Figure 10. The distribution of sini values for stars with M > 0.85 Mq (a) and for stars with M < 0.85 Mq (b). These are compared to 
the expected distribution assuming stars have random rotation axis orientations and that there are no systematic biases in the measured 
parameters (the "No Correction" line), as well as to a model that includes differential rotation (the k = 0.08 line in panel (a)), and 
a model that includes a scale factor a applied to the sini values (the a = 1.24 line in panel (b)). We also show the probability that 
the observed sini distribution is drawn from a model distribution for models invoking a scale factor a (c) and for models that include 
differential rotation (d). 



iTerndrup etafl (2000). A total of 223 of the stars for which 
we detect periods also have a v sin i value given in one of 
these sources. 

Figure[9]compares the photometric rotation periods and 
the v sin i values, and also compares the inferred sin i to the 
stellar mass. The results are shown separately for objects 
selected as lying on the single-star main sequence, and ob- 
jects selected as photometric binaries (see figure 113)) . For the 
remainder of this analysis we only consider nonphotometric 
binaries. The value for sini is determined from: 



i> sin i ■ P 
2nR 



(12) 



where vsini is the measured value of the projected equa- 
torial rotation velocity, P is the measured rotation period, 
and R is the stellar radius. We estimate the stellar radius 
from the Mk magnitude using the 125 Myr ; solar metallicity 
(the Pleiades have [Fe/H]= +0.0 3 ± 0.05; Soderblom et all 



2009), Yonsei-Yale isochrone (Y2; lYi et aT]|200lh . We trans 
form the isochrone to the 2MASS system from the ESQ sys - 
tem using the transformations given bv lCaroenterl (l200lf l. 
To determine the absolute ma gnitudes we assume a dis- 
tance to the Pleiades of 133 pc l|Soderblom et af]|2005l ) and 



an extinction of Ak = 0.01 mag (S07). Note that if we 
instead adopt a distance of 120.2 pc to the cluster base d 
on Hipparcos parallax measurements (|van Leeuwenll2009T ). 
and assume the same extinction and Mx-mass and Mk- 
radius relations, the inferred stellar masses and radii at fixed 
Ks ma gnitude are lower b y ~ 10%. Also note that if we 
use the lBaraffe et all (|l998l ) isochrone models for stars with 
M ~ 0.4 Mq, the inferred masses at fixed Ks are lower by 
~ 8%, while the inferred radii are lower by ~ 30%. The 
inferred radius also depends on the assumed age for the 
cluster. The a s sumed age of 125 Myr for the Pleiades from 
IStauffer et alj (Il998l ) is based on the observed Li-depletion 
of main sequence stars in the cluster, a method which has 
recent ly been called into question for other s tellar popula- 
tions l| Jeffries et al.ll2009l ; lYee fc Jenseiill201Ch . Determina- 
tions of the age of the Pleia des from the main sequence turn- 
off range from ~ 80 My r (lMermilliodlll98ll) to ~ 150 Myr 
IIMazzei fc Piggatdli"989T ). Adopting a younger age of 80 Myr 
for the cluster will increase the radii of stars at fixed mass; 
the affect on the Mx-radius relation, however, is less sig- 
nificant, as stars with fixed mass will also be brighter in 
Mk at younger ages. By comparing the Mx-radius relation 
obtained from the 80 Myr solar-metallicity isochrone from 
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Figure 9. Top: The spectroscopically determined vsini taken 
from the literature vs. the photometric rotation period deter- 
mined in this paper for 223 stars. Open circles show objects 
selected as photometric binaries and grey-scale Xs show non- 
binaries. Arrows show stars for which only upper or lower lim- 
its on v sin i are available. Bottom: sini calculated from vsini, 
the rotation period, and the radii inferred from an isochrone, is 
plotted against stellar mass inferred from an isochrone (see sec- 
tion [32] for details). Stars falling above the line sini = 1 have 
measured rotation period, v sin i and stellar radius values that are 
not fully consistent with one another. The error-bars include un- 
certainties in usini, the period, and the radius (propagated from 
the uncertainty on K s ). We assume a 10% uncertainty on vsini 
for stars without vsini uncertainties given in the literature. We 
exclude known and suspected spectroscopic binary systems from 
both plots. 



ISiess et alj (|2000h using the iKenvon fc Hartmannl (|l995l ) 
temperature-colour relations, to the 125 Myr isochrone from 
the same modelfl, we estimate that the systematic error in 
the radius that may result from overestimating the age of 
the cluster is < 2% over the magnitude range of interest. 
Finally if we use the Mv magnitude to determine the radii 
from the Y2 isochrones, rather than the Mk magnitude, the 
inferred radii for most stars are lower by up to ~ 25%. 

As seen in Figure [9] the sin i values determined with 



We use the ISiess et al.l i2000h isochrones rather than the Y2 
isochrones because the Y2 isochrones do not cover the pre-main 
sequence stage of stellar evolution. 
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Figure 11. a: Photometric rotation frequency (1/P) vs. the fre- 
quency inferred from v sini. The solid diagonal line shows the up- 
per limit (sini = 1). The stars are divided into three mass bins. 
Low frequency (long period) stars with M < 0.85 Mq are more 
likely to lie above the upper limit than higher frequency stars 
with comparable masses. Stars with M > 0.85 Mq generally fall 
below the upper limit for all frequencies, b: Same as above, here 
the velocities of stars with M < 0.85 Mq have been corrected 
with A(To = 0.16kms _1 , and the radii of these stars have been 
increased by 10%, see the text for details, c: The probability that 
stars with M < 0.85 Mq have the same sini distribution as stars 
with M > 0.85 Mq as a function of the applied Aco correction. 
The solid line is for the case when no radius correction is applied, 
the dashed line shows the case when a 10% correction is applied 
to the radii of stars with M < 0.85 Mq. d: Same as C, here a 
A'usini correction is applied. 



equation [12] are greater than 1 for many stars with M < 
0.85 Mq . In Figure [10] we compare the observed distribu- 
tion of sini values to various model distributions. To cal- 
culate the model distributions we generate a sample of sin i 
values assuming the rotation axes of the stars are randomly 
oriented in space (cos i is uniformly distributed) . We then as- 
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sign to each sample sin i value a gaussian error A sin i taken 
from the set of errors associated with the actual data, and 
use these to generate a simulated sample of measured sin 1m 
values with the condition sini m > 0. As discussed in sec- 
tion 13.21 differential rotation may cause the measured rota- 
tion periods to be systematically longer than the equatorial 
rotation periods. This effect will in turn cause the measured 
sin i m values to be systematically greater than the true val- 
ues. To include this effect in our models we associate with 
each simulated sin i value a latitude f3 for the dominant spot 
group. We assume that these spots may be uniformly dis- 
tributed over the surfaces of the stars for reasons outlined 
in section 13.21 so we draw the values of f3 from a uniform 
distribution in sin/3. We then scale each sini m value by the 
factor 1/(1 — fcsin 2 /3) (equation [5]). Finally we parameter- 
ize any remaining systematic errors in the measured values 
of sin i by multiplying the sin i m values by a constant scale 
factor a. We compare the models to the observations us- 
ing the Kolmogorov-Smirnov test, recording the probability 
that the observed sample is drawn from the same distribu- 
tion as the simulated sample as a function of the parameters 
k and a (figure [TOl panels (c) and (d)). 

For stars with M > 0.85 M©, we find that the k = 0, 
q = 1 distribution is consistent with the observed distri- 
bution (the K-S probability that the observed and model 
data-sets are drawn from the same distribution is ~ 50%), 
while a k = 0.08, a = 1 model (i.e. a model with differential 
rotation that is comparable to what is expected for young 
rapidly rotating stars) provides a slightly better match (90% 
probability). A similar result has recently been found in 
an independent study of the Pleiades and Alpha Persei by 
I Jackson &: Jeffries] (|2010l ) who placed limits on the degree to 
which the rotation axes of stars in these clusters are aligned. 

For stars with M < 0.85 Mq the probability that the 
k = 0, a = 1 model is drawn from the same distribution as 
the observations is only ~ 10 -13 . In this case a significant 
value of a = 1.24 ±0.05 is required to fit the observations. If 
we fix a = 1.0 and allow k to vary, the best matching model 
has a very high value of k = 0.60, with a probability of being 
drawn from the same distribution as the observations of less 
than 30%. 

While the distribution of sini values for M > 0.85 Mq 
stars are consistent with the stars having random orien- 
tations and perhaps having a slight degree of differential 
rotation, the measured values of sini for stars with M < 
0.85 Mq appear to be systematically too large (the mode 
is at sini > 1). The combination of measured parameters 
Pv sini/7? may be systematically larger by a factor of ~ 1.24 
than the same combination using the real physical parame- 
ters would be. This suggests systematic errors in one or more 
of the parameters. Below we consider each of the parameters 
in turn. 

Radius - While it is well-known that the radii of rapidly 
rotating stars with M < 0.85 ar e systematically larger than 
theoretical models predict (e.g. I Fernandez et al ] 120091 and 
references therein), the discrepancy appears to be ~ 10% 
and not 24%. If we a dopt the Hipparcos d istance of 120.2 pc 
to the cluster, use the lBaraffe et all (1 19981 ) isochrones rather 
than the Y2 isochrones, or use the V magnitude rather than 
the K magnitude to determine the radii, the inferred radii 
of the stars would be even smaller, exacerbating the sini 
problem. 



Period - While an extreme differential rotation law of 
k = 0.60 could fit the distribution, the match is still not very 
good. Alternatively if the dominant spot groups on lower 
mass Pleiades stars are not randomly distributed, but are 
instead preferentially located at high latitudes, this would 
cause the measured sini values to be systematically larger 
than they would be if the spots are randomly distributed. 

v sin i - The v sin i measurements of cool Pleiades dwarfs 
could be systematically biased toward larger values. Fig- 
ure lllb compares the rotation frequency determined from 
the photometric period to the rotation frequency inferred 
from vsini (i.e. v smi/2irR). For stars with M < 0.85 Mq, 
it appears that lower frequency (longer period) stars are 
more likely to have sin i > 1 than higher frequency (shorter 
period) stars. An offset error in v sin i would yield this effect. 
Most of the longer period, low-m ass stars with sini > 1 have 
v sini measurements taken from lQueloz et all (1 19981 ). These 
authors calculate v sin i from the broadening of a spectral 
cross-correlation function using the relation: 

v sin i = Ay/ a 2 - a\ (13) 

where a 2 is the measured width of the cross-correlation func- 
tion and A and a 2 are parameters which are calibrated by ar- 
tificially broadening the observed spectra of slowly rotating 
field stars to a fixed value of v sin i, and measuring the widths 
of the resultin g cross-correlation fu nctions. For the ELODIE 
spectrograph. iQueloz et all (|l998T ) give A = 1.9 ± 0.1, and 

a = 0.27(B - V) 2 ± 4.51(±0.06) (14) 

where B — V is the measured color of a star, and ao is 
measured in kms . A systematic error in ao would impact 
slower rotators more significantly than faster rotators. To 
estimate the required systematic error in ao, we adjust the 
velocities of stars with M < 0.85 Mq by 

. . A 2 a Aa , . 

Aosin! = (15) 

v sin i 

and use the K-S test to compare the resulting sin i distribu- 
tion to the sini distribution for stars with M > 0.85 Mq. 
For this test we assume <ro = 4.8 for these stars. Figure [TTt 
shows the probability that the two samples are drawn from 
the same distribution as a function of Aao. We conducted 
a similar test assuming a constant Aw sini (figure [TTJI) . We 
find Aa = 0.6±0.2kms _1 or Awsini = -1.5 ± 0.5 kms -1 
yields a sin i distribution for stars with M < 0.85 A/0 that is 
statistically indistinguishable from the distribution for stars 
with M > 0.85 Mq. If the isochrone radii of stars with 
M < 0.85 Mq are also assumed to be 10% too low, we 
find Aao = 0.33lo;2?kms _1 or Ausini = -0.85 ±0.63. Fig- 
ure lllb shows how a correction of Aao = 0.33, together with 
a 10% radius correction affects the frequencies inferred from 
v sin i for stars with M < 0.85 Mq. A systematic error of 
Aao = 0.33^q' 27 kms -1 is consistent at the la level with the 
systematic uncerta inty in ao of ~ 0.06 kms -1 estimated by 
IQueloz et al.l |l998). Differences between the young Pleiades 
stars and the older field stars used to calibrate the (B — V)- 
ao relation could result in an even grea ter systematic error in 
ao. For example. IStauffer et alj (|2003l ) showed that Pleiades 
stars do not have typical spectral energy distributions-in 
particular they have excess emission in the B band (and 
are therefore bluer in B — V than expected from their ef- 
fective temperatures) . Equation [14] may therefore underes- 
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timate the intrinsic broadening of Pleiades K dwarfs. Quan- 
titatively, a 0.1 mag shift in B — V for Pleiades stars with 
B — V ~ 1.0 yields an additional ~ 0.05 kms -1 shift in co- 
in conclusion, assuming that the radii are systemat- 
ically underestimated by ~ 10%, and that the intrinsic 
broadening of these stars are also underestimated by ~ 
0.1 kms -1 , we find that the observed distribution of sini 
values would then be consistent with the model distribution 
for a differential rotation parameter of k = 0.27 ± 0.14 (the 
M < 0.85 M@, a = 1.1, Aa = 0.1 line in figure [TOhV Al- 
lowing for a larger error in the intrinsic broadening would 
reduce the required k. If we only apply the radius correc- 
tion, the differential rotation parameter would need to be 
rather large (k = 0.38±g:^; the M < 0.85 Mq, a = 1.1 line 
in figure HHt O. We conclude that some combination of the 
above effects provides a plausible explanation for the large 
number of low-mass stars with measured sini > 1. 

Finally, we note that a qualita ti vely s imilar effect has 
also been noted by IJackson et al.l (|2009h for stars with 
0.2 Mq < M < 0.7 Mq in the comparably aged open 
cluster NGC 2516. They find that the radii of stars with 
M ~ 0.2 Mq must be ~ 50% larger than theoretical predic- 
tions for the v sin i and period data to be consistent in this 
cluster. 



4.3 The Period-Mass Relation and Comparison to 
Other Clusters 

Figure [12] shows the rotation period vs. absolute magnitude 
Mk and mass, together with the relation between v sin i 
and absolute magnitude Mk /mass. By comparing the v sin i 
values of stars without photometric period detections to the 
?;sini values of stars with period detections, it is apparent 
that over the mass range 0.6 Mq < M < 1.0 Mq there 
does not appear to be a bias against detecting photometric 
periods for slow rotators. For star with M > 1.0 Mq there 
does appear to be a bias against detecting periods for the 
slowest rotators, so our sample is not complete in this mass 
range. For stars with M < 0.6 Mq the bias appears to be 
against detecting rapid rotators. This is likely due to short 
period stars being redder on average than long period stars 
(see Fig. 1131 which demonstrates this for (V — Ks), this is 
likely to be true for (r — Ks) as well), so that at fixed Ks 
short period stars are fainter in r than long period stars, 
and therefore have poorer precision HATNet light curves. 

In figure [T3] we compare the period-mass relation 
for the Pleiades to the relations for four other simi- 
larly aged clusters. These include 3 clusters studied by 
the M ONITOR projec t (NGC 2547, 40 Myr, llrwin et all 
|200S| ; M 50, 130 Myr, llrwin et all 120091 ; and NGC 2516, 
150 Mvr.lfrwin et al.ll2007r i. and the cluster M35 (180 Myr, 
iMeibom et"al]|2009h . For the MONITOR clusters we adopt 
the stellar masses given in their tables. For M35 we de- 
termine stellar mas ses using the extinc tion corrected V- 
magnitudes given bv lMeibom et all (|2009T ) tog ether with the 
Y2 is ochrones. We assum e an age of 180 Myr l|Kalirai et al.l 
2003), distance of 912 pc (Kalirai et al. 2003), and metallic- 
ity of [Fe/H] = -0.21 (ISchuler et al.ll2003h . 

The data presented here can be combined with data 
for other open clusters to test theories of stellar angular 
momentum evolut ion. While a soph i sticate d analysis like 
that presented by iDenissenkov et all ((2009) is beyond the 
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Figure 13. (V — Ks) vs. Mk CMD for probable Pleiades mem- 
bers with detected rotation periods. The color-intensity of each 
point is scaled by the star's rotation period. At fixed Mjf, short 
period rotators tend to be redder in (V — Ks) than long period 
rotators. This effect was previously noticed using v sin i data by 
IStauffer et"afl ||2003|) for Pleiades stars an d bv lAn et all h00l\) for 
stars in NGC 2516. IStauffer etaD HoojJ suggested that this ef- 
fect may be due to rapid rotators having more red spots than slow 
rotators. Alternatively this may be due to strong magnetic fields 
in rapid rotators inhibiting convection, which ca uses these stars 
to be larger and cooler than slow rotators (e.g. [Chabrier et al] 
2007). Stars located above the solid line are selected as poten- 
tial photometric binaries, and are excluded from the analysis in 
section 14.21 



scope of this paper, it is instructive to compare the Pleiades 
sample to M35. As seen in figure [T5l the main period-mass 
sequences in the two clusters overlap, indicating that slow 
rotators with 0.7 Mq < M < 1.1 Mq do not spin-down 
between 125 Myr and 180 Myr. We also show in figure [15] 
the substantial expected evolution between the two clusters 
for a standard soli d-body rotation angul ar-momentum evo- 
lution model (e.g. Irlartman et al.|[2009al ). While the slowly 
rotating stars do not spin-down between the ages spanned by 
the two clusters, more rapid rotators with M > 0.7 Mq do 
appear to have spun-down. These two features can be repro- 
duced by models that invoke core-envelope decoupling with 
a coupling time-scale that depends on the rotation period, 
such that rapid rotators have a short coupling time-scale 
of a few Myr, while slow ro t ators have a long time- scale 
of ~ 100 Myr (iBouvieri 120081 ; IDenissenkov et al.|[2009T ). As 
noted by iBouvi er (2008), one consequence of slow rotators 
having a less efficient core-envelope coupling than rapid ro- 
tators is that slower rotators should exhibit more significant 
Li depl etion than rapid rotato rs. Such an effect has been 
soon bv lSoderblom et all l|l993h for the Pleiades. 



5 SUMMARY 

In this paper we have presented new rotation period mea- 
surements for 368 stars in the Pleiades, increasing by a fac- 
tor of 5 the number of stars in this cluster with measured 
periods. This includes 14 newly identified probable cluster 
members. The sample is ~ 93% complete for stars with 
0.7 Mq < M < 1.0 Mq. By comparing our sample to the 
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Figure 12. Left: photometric rotation period vs. the Kg absolute magnitude/mass for 350 probable members of the Pleiades cluster. 
Note that the period is plotted on a logarithmic scale. Filled circles show the adopted period values. Crosses show the periods prior 
to applying the alias correction described in section 13.11 Open circles show the periods after this correction. Right: spectroscopically 
measured vsiai vs. the Kg absolute magnitude/mass for stars in our sample with measured photometric rotation periods (filled circles) 
and stars without photometric rotation periods (crosses). Note that risini is plotted on an inverted logarithmic scale. For stars with 
0.6 Mq < M < 1.0 Mq the sample with measured periods is nearly complete, and there does not appear to be any significant bias 
toward selecting fast rotators. For stars with M > 1.0 Mq there does appear to be a bias toward detecting photometric periods for 
faster rotators. For stars with M < 0.6 Mq there may be a bias against selecting shorter period rotators. This is likely due to rapid 
rotators having slightly redder (V — Kg) colors at fixed Kg magnitude (fig. 1131 1. making fast rotators slightly fainter in the r filter used 
by HATNet than slower rotators at fixed Kg. 



large sample of Pleiades stars with v sin i measurements in 
the literature, we find that our sample is not biased toward 
short periods over this same mass range. 

We show that for stars with M > 0.85 Mq the inferred 
distribution of sini values is consistent with the stars hav- 
ing an isotropic distribution of rotation axes at the la level, 
the agreement between the model and observations is even 
better if a k — 0.08 differential rotation law is assumed; a 
differential law of this form is consistent with theoretical 
expectations for zero-age main s equence stars. Thi s resul t 
is consistent with the findings of I Jackson fe Jeffries! l|2010h . 
For stars with M < 0.85 Mq the inferred sin i values are sys- 
tematically larger th an 1, and is similar to what was seen by 
IJackson et al.l (2009) for stars in NGC 2516. Our observa- 
tions imply that the combination P(v sini)/ R is too large by 
~ 24%. We argue that a ~ 10% systematic error in the radii 
of these stars, together with a > 0.1 kms -1 systematic error 
in the assumed intrinsic broadening of their spectral lines, 
and a k ~ 0.3 differential rotation law provides a plausible 
explanation for this discrepency. 

Finally, we find that the mass-period d iagram for the 
Pleiad es is remarkably similar to that seen bv lMeibom et all 
(2009) for the ~ 185 Myr cluster M35. In particular there 
appears to be very little change in the rotation rates of the 
slowest rotators with 0.7 Mq < M < 1.1 Mq between the 
ages of these clusters. This result provides support for claims 
that these stars have inefficient internal angular momentum 
transport, and exhibit significant core-envelope decoupling. 
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Table 1. Catalogue of rotation periods for Pleiades stars. The full table includes 383 stars - 368 have periods determined from our photometry, and 15 have periods taken from the 
literature. We exclude the star HII 3167 with a period listed on WEBDA because it is not a cluster member, and we include the star CFHT-PL-8 which is not included in the S07 
catalogue. We also include 4 stars excluded from the S07 catalogue, which have previously been identified as candidate cluster members, and which we found to be probable cluster 
members in section [3.51 and 14 stars which we have newly determined to be probable cluster members. The former stars are SRS 34337, SRS 79807, HCG 84 and HCG 235, while 
the latter stars can be identified by having no "OtherlD" given. The first ID is an internal HAT- ID and is the ID used for labeling the light curves, the other IDs are defined in S07. 
The formal error on the period from fitting a Gaussian to the peak in the periodogram is given by ePer. The lcr width of the Gaussian peak (in frequency) is given by fwidth. Other 
references with periods for the star in question are given in rPer, the actual periods are given in Table [4] The r-band peak-to-peak amplitude of the light curve is given by rAmp; stars 
without an amplitude listed have period s take n from the literature. T he mass and radius ar e estimated from the K magnitude and the YY2 isochrone s as described i n the text, except 
for 9 faint stars from IScholz fc Eis loffel (2004) and 1 faint star from iTerndrup et al . (1999) for which the masses and radii are determined using the iBaraffe et al isochrones. 
The reference for the vsini value is given by rvsini, and fbinary is a flag indicating if the star is a known spectroscopic or visual binary system, it also includes the reference for this 
determination. Positions and magnitudes are taken from S0 7. Stars which were not included in the S07 catalogue have positions and magnitudes take n from 2MASS or the 2MASS 6x 
PSC. For CFHT-PL 8 we take the optical photometry from lSamus et al.l d2003T ) , for SRS 34337 and SRS 79807 we take the optical photometry from lSchilbach et al.l lll995l) . while for 
HAT-259-0000829 we take the optical photometry from the Tycho-2 catalogue (|Hog et al.ll200Q|y The full table is included in the supplementary material of the online edition of the 
journal. A portion is shown here for guidance regarding its form and content. Abreviations for the references are defined in the ReadMe file associated with the electronic version of the 
table. 



ID 


RAdeg 


DEdog 


OtherlD 


Nobs 


RMS 


Bmag 


Vmag 


Kmag 


Mass 


Radius 


vsini 


rvsini 


fbinary 




J2000 


J2000 






mag 


mag 


mag 


mag 


M 


Rq 


km s — 1 






HAT-259-0000604 


51.925262 


23.803688 


PELS121 


6047 


0.01314 


10.96 


10.30 


8.679 


1.150 


1.057 


4.900000 


Q98 




HAT-259-0001639 


52.168613 


25.607782 


AKIII59 


6020 


0.00581 


12.59 


11.75 


9.723 


0.832 


0.732 


-99.990000 






HAT-259-0002068 


52.409874 


24.510546 


DH008 


6055 


0.02236 


-9.99 


-9.99 


9.698 


0.839 


0.738 


-99.990000 






HAT-259-0001310 


53.274277 


22.134232 


PELS004,DH026,AKII131 


6046 


0.00632 


12.25 


11.43 


9.425 


0.911 


0.809 


2.600000 


Q98 




HAT-259- 0010017 


53.293736 


22.522045 


DH027 


5983 


0.04235 


-9.99 


-9.99 


11.248 


0.526 


0.502 


-99.990000 






HAT-259-0000347 


53.530495 


24.344501 


PELS006,TrS4,AKIII288 


6035 


0.01747 


10.08 


9.57 


8.274 


1.309 


1.247 


35.900000 


Q98 




HAT-259-0004453 


53.684944 


26.096003 


DH037 


5977 


0.03269 


-9.99 


-9.99 


10.371 


0.683 


0.597 


-99.990000 






HAT-259-0000705 


53.696770 


26.094885 


PELS007,AKIII327 


6009 


0.00490 


11.14 


10.48 


8.833 


1.097 


1.000 


2.700000 


Q98 




HAT-259-0000440 


53.882030 


22.823627 


PELS124 


6050 


0.01401 


10.40 


9.86 


8.541 


1.200 


1.115 


20.500000 


Q98 




HAT-259- 0001134 


54.004112 


24.266083 


AKIII391 


6013 


0.00495 


12.06 


11.29 


9.388 


0.922 


0.820 


2.300000 


Q98 





Table 2. Catalogue of Pleiades members observed by HATNet without period determinations. The full table includes 120 stars. We include the stars HII 708, HII 727, HII 975, and 
HHJ 409, which have period measurements taken from the literature, and are included in table 1, but did not have periods recovered from the HATNet light curves. The columns 
are as in table 1, here we also include Nobs, which is the number of photometric observations, and the unbiased RMS of the HATNet TFA light curve for the target. Only stars with 
9.5 < r < 14.5 are included in this table. The full table is included in the supplementary material of the online edition of the journal. A portion is shown here for guidance regarding its 
form and content. Abreviations for the references are defined in the ReadMe file associated with the electronic version of the table. 





ID 


RAdeg 


DEdog 


Per 


rAmp 


,7 mag 


Hmag 


Kmag 








J2000 


J 2000 


d 


mag 


mag 


mag 


mag 


HAT- 


•213- 


-0004432 


46.376495 


27.295223 


1.47490588 


0.0421 


11.386 


11.002 


10.921 


HAT- 


•213- 


■0000301 


46.395025 


26.326672 


0.49278549 


0.0068 


7.894 


7.317 


7.198 


HAT- 


•213- 


-0003310 


46.423947 


26.964050 


0.51638485 


0.0301 


10.947 


10.575 


10.448 


HAT- 


•258- 


■0007599 


46.452067 


26.214544 


0.80443471 


0.0680 


11.960 


11.292 


11.068 


HAT- 


•258- 


■0001512 


46.545842 


24.555460 


1.56383216 


0.0114 


10.440 


10.123 


10.063 


HAT- 


•213- 


-0000501 


46.549253 


26.699934 


1.02287839 


0.0544 


7.864 


6.928 


6.573 


HAT- 


•258- 


-0005170 


46.567769 


24.944794 


7.96319798 


0.0110 


11.796 


11.298 


11.198 


HAT- 


•213- 


-0005056 


46.574536 


26.634954 


6.48431835 


0.0192 


11.557 


11.194 


11.097 


HAT- 


•258- 


-0002836 


46.593644 


24.553173 


0.12926165 


0.0015 


11.035 


10.531 


10.453 


HAT- 


•213- 


-0008723 


46.597723 


26.747597 


0.61090482 


0.0228 


12.609 


12.438 


12.416 



Table 3. Catalogue of periodic variables in HATNet field G259 that are not cluster members. The full table includes 1804 stars. The positions and magnitudes are taken from 2MASS. 
The amplitude is the peak-to-peak amplitude of the best-fit sinusoid to the signal-recovery-mode TFA light curve of each star. Note that the true amplitude may be systematically 
higher by a factor of ~ 1.3 (see section r3.3l l. The full table is included in the supplementary material of the online edition of the journal. A portion is shown here for guidance regarding 
its form and content. 



ID 


Period 

d 


Rof 


BPL106 


0.170000 


S04 


BPL115.DH476 


0.121250 


S04 


BPL125,DH494 


0.806250 


S04 


BPL138.DH505 


1.075420 


S04 


BPL150 


0.769167 


S04 


BPL164.DH544 


0.840000 


S04 


CFHT-PL8 


0.401000 


T71 


HCG20,T1B,DH105 


2.700000 


K98 


HCG254,BPL129 


0.401667 


S04 


HCG346,HHJ111,BPL190 


1.677920 


S04 



Table 4. Previous rotation period measurements for Pleiades stars. The full table includes 94 measurements for 66 stars. The full table is included in the supplementary material of the 
online edition of the journal. A portion is shown here for guidance regarding its form and content. The references are defined in the ReadMe file associated with the electronic version 
of the tabic. 



